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COLLIDING ELECTROMAGNETIC WAVES: 
THE BELL-SZEKERES SOLUTION 

In an important paper, Bell and Szekeres (1974) gave an exact solution 
which describes the collision and subsequent interaction of two electro- 
magnetic plane waves. This solution, which appears to be remarkably 
simple, will be described in this chapter. Other exact solutions describing 
colliding electromagnetic waves will be described in Chapters 16 and 17, 
after more powerful techniques have been developed. 

15.1 The Bell— Szekeres solution 

Bell and Szekeres have considered a very simple situation involving a 
collision of two step electromagnetic waves. Prior to the collision, the 
waves may be described by the familiar line element (4.15). In regions I 
and II, this takes the form 

ds 2 = 2dwdr + a 2 Q(u)(X 2 + Y 2 )du 2 - dX 2 - dY 2 (15.1) 

and the step wave is given by $22 = a?Q(u). The opposing wave with 
^oo = b 2 Q(v ) may be described in an identical way with the null coordi- 
nate u replaced by v, and a different space-like coordinate replacing r. 

The two electromagnetic field components in the initial regions I, II 
and III are taken to be 

$ 2 =a0(w), and $ = bO(v). (15.2) 

The fact that both of these expressions are real simultaneously indicates 
that the polarization of the two electromagnetic waves is aligned. In 
addition, they have identical step profiles. The relative amplitudes a 
and b may, of course, both be equated to unity when convenient. The 
components of the Weyl tensor are initially all zero. 

By this stage, we are familiar with the fact that it is appropriate to 
transform the metrics in all regions into Rosen form. Accordingly, the 
line elements in the initial regions may be taken in the forms of (4.19) 
and (4.20), namely 

I (u < 0, v < 0) ds 2 = 2dudv - dx 2 - dy 2 

II(u>0,v<0) ds 2 = 2dudv - cos 2 au(dx 2 + dy 2 ) (15.3) 

III (u < 0, v > 0) ds 2 = 2dudv - cos 2 bv(dx 2 + dy 2 ). 
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The initial- value problem is now well set, and it remains to find the unique 
solution in the interaction region. 

The field equations for colliding electromagnetic waves have already 
been obtained in Chapter 6. In terms of the metric functions of the 
Szekeres line element (6.20), the gravitational field equations for the in- 
teraction region take the form (6.22), and Maxwell's equations are given 
by (6.21). The appropriate boundary conditions for this situation have 
been described in Chapter 7. 

As usual, it may be noted that (6.22a) can immediately be integrated 
to give 

e- u = f(u)+g(v) (15.4) 

where / and g are arbitrary decreasing functions in the interaction re- 
gion. From the junction conditions it can be seen that, in this case, these 
functions necessarily take the forms: 

f = \ — sin 2 au, g = | — sin 2 bv. (15.5) 

The solution of the complete set of equations (6.21,22) given by Bell and 
Szekeres is 

U = — log cos(au — bv) — log cos(au + bv) 
V = log cos(au — bv) — log cos(au + bv) (15.6) 
W = 0, M = 0, $ 2 = a, $ = b. 

It can be seen that these functions satisfy the required O'Brien-Synge 
boundary conditions. Thus the line element in the interaction region IV 
is 

ds 2 = 2dudv — cos 2 (au — bv)dx 2 — cos 2 (au + bv)dy 2 . (15-7) 

This line element is in fact one form of the Bertotti-Robinson solution 
(Bertotti 1959, Robinson 1959), which is known to be conformally flat. 
The coordinate transformation 

t + r = coth ( -sech -1 cos(att + bv) — — | 

t — r = — tanh f -sech -1 cos(cm + bv ) + — ) (15.8) 

7T X 

9 =au — bv H — , d> = — 

2 Y 2 

transforms the line element (15.7) into the more familiar form 

a 2 

ds 2 = ^(dt 2 - dr 2 - r 2 d6 2 - r 2 sin 2 Odcf 2 ). (15.9) 
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Halilsoy (1987) has considered an application of the transformation 
u — > u'(u), v — > v'iy) to the above solution. However, this does not lead 
to any new results that are of physical significance. 

15.2 The structure of the solution 

The above solution describes the collision of two step electromagnetic 
waves whose polarization vectors are aligned. This follows from the fact 
that $02 is real, and W = 0. 

As has been pointed out above, the solution inside the interaction 
region IV is conformally flat. However, as Bell and Szekeres have shown, 
there are necessarily discontinuities in the derivatives of the metric func- 
tion V across the initial boundaries of this region. These manifest them- 
selves as the impulsive gravitational waves 

$ 4 = - a tan bvS(u)Q(v), ^ = -b tan au 5 (v)G(u) (15.10) 

which may be considered to be generated by the collision. 

It is in fact a general feature of colliding electromagnetic plane waves 
that gravitational waves are always generated by the collision. This can 
easily be demonstrated by considering the field equations (6.22d,e). In the 
interaction region, $2 and $0 are necessarily both non-zero. It follows 
that V and W can not both remain constant. Thus, the congruences 
tangent to both waves must begin to shear, and it can be seen from (6.23) 
that the components of the Weyl tensor * and $4, and possibly $2, will 
necessarily appear. In this case, the discontinuity in the electromagnetic 
components causes impulsive gravitational waves to be generated. It can 
also be shown that smooth-fronted electromagnetic waves would generate 
smooth-fronted or step gravitational waves that wouls persist through the 
interaction region. 

It may also be observed from (15.10) that difficulties occur at the 
points u = 11 /2a, v = and v = it /2b, u = 0. At first sight, it would 
appear that curvature singularities occur at these points (Matzner and 
Tipler, 1984). However, curvature tensors of this type can be interpreted 
as distributions (see Geroch and Traschen, 1987) and, in this case, the 
standard definition of a curvature singularity in terms of unboundedness 
in a parallely propagated frame is inapplicable. 

The global structure of this solution has been analysed in detail by 
Clarke and Hayward (1989). They have confirmed that, in regions II and 
III, the surfaces on which / = —1/2 and g = —1/2 respectively behave as 
'fold' singularities similar to those of the Khan-Penrose solution described 
in Section 8.2. 
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Initially, the most surprising feature of the Bell-Szekeres solution 
was that a space-like curvature singularity does not occur in region IV. 
In this case, the surface on which / + g = 0, which is here given by 
au + bv = it/2, turns out to be merely a coordinate singularity since the 
curvature tensor is clearly bounded. Bell and Szekeres have demonstrated 
how this singularity can be removed by considering the following coordi- 
nate transformation: 

T = — cos(au + bv) cosh cy 

Z = cos(au + bv) sinhcy 

1 ' (15.11) 

X = cos(au — bv) cos cx 

Y = cos(au — bv) smcx 
where c = \ / 2ab. With this, the line element (15.6) becomes 
2 / (1 + Z 2 )df 2 - 2TZdfdZ - (1 - f 2 )dZ 2 ' 



ds' 



c 2 (l —f 2 + Z 2 ) 

(1 - Y 2 )dX 2 + 2XYdXdY + (1 - X 2 )dY 2 
c 2 (l -X 2 -Y 2 ) 



(15.12) 



which is clearly regular when T = Z = 0. 

In order to retain the one-to-one correspondence, it is appropriate ini- 
tially to restrict the Bell-Szekeres coordinate x to the range — n < cx < n. 
The coordinates X and Y are restricted by X 2 +Y 2 < 1, where the bound- 
ary X 2 + Y 2 = 1 contains the space-like plane on which the collision 
occurs. 

It has been shown by Clarke and Hayward (1989) that region IV is 
part of a space-time that is regular apart from two covering-space sin- 
gularities that occur when au — bv = ±7r/2. These are points on the 
boundary of region IV where it joins with the limits of regions II and 
III. At these points, which are on the polar axis 6 = and 9 = ix in the 
Bertotti-Robinson form of the metric (15.9), X = Y = 0. 

With the above transformation, each interval of the coordinate x of 
length 2n/c corresponds to a circular region in the X,Y plane with the 
origin removed and cut along the line Y = 0, X < 0. Neighbouring 
regions of x are joined together along these cuts to form a surface which 
continuously winds around the origin. The origin of the X, Y plane thus 
corresponds to a singularity of a quasiregular covering-space type. 
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It is also sometimes convenient (Clarke and Hayward, 1989) to make 
the further coordinate transformation 




(15.13) 



4> = sin 1 



With this, the line element (15.12) takes the alternative form of the 
Bertotti-Robinson line element 

ds 2 = ^ (cosh 2 p d X 2 - dp 2 - d9 2 - sin 2 9 d^ . (15.14) 

Region IV is now described by the line element (15.14) with the coordinate 
ranges given by 



— OO < (ft < oo 
— oo < p < OO 
— 7r/2 < x < — sin -1 | tanhp| 
cos -1 (cos x cosh p) < 9 < cos -1 (— cos x cosh p). 



(15.15) 



This is joined to regions II and III across the surfaces 

cos 9 = cos x cosh p and cos 9 = — cos x cosh p. (15.16) 

Clarke and Hayward (1989) have also shown that the Bell-Szekeres 
solution is extendable to the future through the surface au + bv = n/2, 
but the extension is not unique because of the presence of the singulari- 
ties. They have further suggested two possible extensions which will be 
described in the next section. 

It is also somewhat remarkable that the global structure of the Bell- 
Szekeres solution, as described by Clarke and Hayward (1989), is very 
similar to that of the collision of an impulsive gravitational wave with a 
null shell of matter, as described by Babala (1987). This will be discussed 
in Section 19.3. 

Chandrasekhar and Xanthopoulos (1988) have conducted a thorough 
perturbation analysis of this particular solution. Throughout the interac- 
tion region subsequent to the collision they have obtained a complete set 
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of bounded normal modes that are expressed in terms of spin-weighted 
spherical harmonics. These modes exhibit a behaviour of ever-increasing 
frequency as the Cauchy horizon is approached. 

Chandrasekhar and Xanthopoulos (1988) have also considered per- 
turbations in the initial regions which contain the approaching waves. For 
the approaching wave in region II which propagates along the null line 
u = const, v < 0, they have shown that ^-independent perturbations are 
not permitted to all orders, while all the independent perturbations ex- 
hibit strong divergences along the 'fold' singularity. Unfortunately, these 
perturbations can not be joined continuously to those of the interaction 
region. It is therefore not possible to consider the effects of these pertur- 
bations subsequent to the collision. 

15.3 Extensions of the solution 

In the Bell-Szekeres solution, it has been shown that the hypersurface 
/ + g = on which the opposing waves mutually focus each other is 
a Cauchy horizon rather than a curvature singularity. It is therefore 
appropriate to consider possible extensions to this solution through the 
horizon, even though any such extension will not be unique. Two possible 
extensions have been given by Clarke and Hayward (1989). 

One natural extension of region IV is simply the future extension of 
the region as described by the line element (15.14). This is achieved by 
extending the coordinate range to 

— oo < 6. p < oo, — 7r/2 < v < oo 

_' „ ' (15.17) 

cos A < 6 < cos (-A) 

where A = min(l, cos x cos p). In this extension, the waves cross and 
focus. The subsequent space-time is unchanging. 

Another equally natural extension of region IV is obtained by noting 
that the line element (15.14) admits a reflection symmetry x — > ~X about 
the hyperplane x = (or T — > — T about T = 0). Thus, there is an 
extension in which region IV is extended as far as x = and is then 
followed by the time reverse of the entire solution. In this extension, 
the waves cross, mutually focus each other, re-expand and then separate 
leaving Minkowski space between them to the future. This situation is 
illustrated in Figure 15.1 in terms of the original null coordinates and the 
line elements (15.3) and (15.7). 

Other apparent extensions of the Bell-Szekeres solution have been 
suggested by Giirses and Halilsoy (1982) in which the approaching elec- 
tromagnetic waves consist of a sequence of steps. The step lengths, how- 
ever, have to be double the distance to the fold singularities. Thus, these 
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Figure 15.1 An extension of the Bell-Szekeres space-time. Regions marked I 
are flat. Regions marked II and III contain electromagnetic waves with $2 = 0,, 
and $0 = b respectively, having metrics given by (15.3). There are impulsive 
gravitational waves along the boundaries of region IV, which has metric (15.7). 

solutions necessarily extend the space-time through the covering space 
singularities and can not, therefore, be considered as appropriate solu- 
tions for colliding plane waves. In all these solutions the space-time can 
only be extended uniquely up to the fold singularities in regions II and 
III and the surface in region IV on which au + bv = 71/ 2. It must there- 
fore be concluded that the Giirses-Halilsoy extensions are not physically 
significant. 

15.4 A non-colinear collision 

In the Bell-Szekeres solution, the polarization of the approaching elec- 
tromagnetic waves is aligned. It is appropriate now to consider how this 
solution may be generalized to include the case when the polarization is 
not aligned. 

In the approaching waves a rotation of the polarization is expressed 
by a rotation of the components $ 2 or <£> i n the complex plane. The 
metrics in both regions are unaltered, since a null electromagnetic field 
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is defined by the metric only up to a constant duality rotation. Thus, to 
extend the Bell-Szekeres solution, the same initial conditions are given 
except that the component $02 becomes complex and, from (6.22 e), W 
must become non-zero in the interaction region. 

A solution presented by Griffiths (1985) satisfying this condition may 
be written in the form 

e~ u = I cos 2au + | cos 2bv 

2V ( 1 + cos 2au cos 2bv + sin 2au sin 2bv cos a) 

(1 + cos 2au cos 2bv — sin 2au sin 2bv cos a) ^ ^ 

, sin 2au sin 2bv sin a 

sinh W = — — 

(cos 2au + cos 2bv) 

M = 0, $ 2 = ae ie , $ = ae _i * 



where 



6 = g cos 1 (cosa sechVF). (15.19) 



The relative polarization angle of the two waves prior to the collision is 
given by a. It may be noticed that the expression for U here is equivalent 
to those of (15.4,5) and (15.6). The common multiple 2 has been inserted 
so that this solution can be more directly related to the Bell-Szekeres 
solution to which it reduces when a = 0. 

It can be shown that this solution has exactly the same properties 
as the aligned Bell-Szekeres solution. It is conformally flat everywhere 
except at the boundaries of the interaction region where there are impul- 
sive gravitational waves that may be considered to be generated by the 
collision. These are given by 

. ' 1 ' (15.20) 
*o = -b tan aue~ ia 5 (v)e(u). 

According to a theorem of Tariq and Tupper (1974), this new solution 
must also be the Bertotti-Robinson space-time in a different coordinate 
system. In fact it can be obtained from the Bell-Szekeres line element 
(15.7) by the simple rotation 

a . a . 
x = cos —x + sin — y 

a , a , ( 15 - 21 ) 

V — — sin —x + cos —y . 
y 2 2 



Thus, the gravitational field in this case is identical to that of the Bell- 
Szekeres solution, as may have been expected. In the collision of plane 
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electromagnetic waves with non-aligned polarization, it must therefore be 
concluded that the dynamics remains unaltered. It is only the form of the 
metric in the interaction region that needs to be modified to that given 
in this section. 



